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ABSTRACT 

We report the discovery of 14 low-mass binary systems containing mid-M to mid-L dwarf 
companions with separations larger than 250 AU. We also report the independent discovery of 
9 other systems with similar characteristics that were recently discovered in other studies. We 
have identified these systems by searching for common proper motion sources in the vicinity of 
known high proper motion stars, based on a cross-correlation of wide area near-infrared surveys 
(2MASS, SDSS, and SIMP). An astrometric follow-up, for common proper motion confirmation, 
was made with SIMON and/or CPAPIR at the OMM 1.6 m and CTIO 1.5 m telescopes for all 
the candidates identified. A spectroscopic follow-up was also made with CMOS or GNIRS at 
Gemini to determine the spectral types of 11 of our newly identified companions and 10 of our 
primaries. Statistical arguments are provided to show that all of the systems we report here 
are very likely to be physical binaries. One of the new systems reported features a brown dwarf 
companion: LSPM J1259-I-1001 (M5) has an L4.5 (2M1259-I-1001) companion at ^340 AU. This 
brown dwarf was previously unknown. Seven other systems have a companion of spectral type 
LO-Ll at a separation in the 250-7500 AU range. Our sample includes 14 systems with a mass 
ratio below 0.3. 


Subject headings: binaries: general — stars: low-mass, brown dwarfs 


1. INTRODUCTION 

In the last two decades, the search for low- 
mass stars and brown dwarfs has intensified with 
the advent of wid e area su rveys such as 2MAS S 
(iGutri et al.ll200.‘^. DE NIS (lEpchtein et al.lll999ll. 
SDSS ( Ahn et al. 2012ll . UKIDSS ( Lawrence et al 


2007), and more recently WISE (|Cutri fc et al 
20121) . In addition, a handful of surveys were 
aimed at finding nearby binary systems in- 
cluding low mass stars such as the sur veys of 
Reid fc GizisI ( 19971) and I Reid et al.l ( 2001 ) for ex- 
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ample. IPhital et al. ( 2010h assembled the SLoW- 
PoKES catalog which contains 1342 binary sys¬ 
tems with a projected separation of > 500 AU, in 
which one of the components is a mid-K to mid- 
M dwarf. Numerous wide binary systems have 
also been discovered individually. Of all the sys¬ 
tems known, very few have a separation larger 
than 7000 AU. The system composed of NLTT 
20346 (an M5V and M6V binary) and 2MASS 
J0850359-I-105716 (an L6 dwarf) i s one of these, 


with a separation of 7700 AU (IRadigan et al 


2009), while Konigstuhl 3 A-BG (IFahertv et al 


20111 ). a hierarchical system composed of an E8V 


orbited by an M8V and L3 binary, has a separa- 
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tion of 12 000 AU. 


An interest of binary systems is that both com¬ 
ponents share the same age and metallicity. For 
systems comprising an M dwarf companion for ex¬ 
ample, this can hel p calibrate the meta llicity scale 
for this type of star ( Bonfils et al.ll2005[) . The com¬ 
ponents of wide visual binaries can also be studied 
in detail individually using seeing-limited instru¬ 
ments, making their characterization much easier. 
Additionally, the evolution of each component of 
a wide binary system has not been influenced by 
the other component, and both can be viewed as 
an “isolated” star. 


The discovery of numerous low-mass stars and 
brown dwarfs has led to a growing interest in un¬ 
derstanding the formation mechanism of these ob¬ 
jects whose masses are b elow the typical Jeans 
mass in a molecular cloud ( Mever et al.ll200r)ll , and 
numerous hypotheses about their origin have been 
formulated over the years. Binary systems con¬ 
taining a low mass component provide a good tool 
for attempting to discriminate between the various 
possibilities, as their properties (e.g., separations, 
mass ratios) ma y be affected differently by dif¬ 
ferent processes ( Duchene k. Krau3l2013ll . Wide 
binary systems with separation larger than hun¬ 
dreds of AU are of particular interest for this pur¬ 
pose as, being more weakly gravitationally bound 
and spanning sizes larger than typical protostel- 
lar disks, they bring additional and perhaps more 
stringent constraints on the underlying formation 
processes. 


In this article, we report the discovery of 14 new 
binary systems that all have a component with 
a spectral type later than M6 and a separation 
larger than 250 AU. We also recovered nine previ¬ 
ously known binaries with similar characteristics. 
Section [5] details our archival search for wide bina¬ 
ries and the identification of our candidates, while 
Section [3] describes the follow-up photometry and 
the calculations of the proper motion. The spec¬ 
troscopic follow-up as well as the results extracted 
from these data are explained in Section |4l The 
probability of random alignement has been calcu¬ 
lated for each system and is presented in SectionjS] 
A short discussion is presented in section [5] Our 
results are summarized in section [71 


2. SEARCH AND IDENTIFICATION 
OF THE CANDIDATE COMPAN¬ 
IONS 


2.1. Sample 


We searched for new binary systems comprising 
a low-mass component through common proper 
motion based on multi-epoch wide area near- 
infrared surveys. As a starting point for po- 
tential primary stars , we consider ed the NLTT 
( Salim fc Gould! 2003) and LSPM ( Leninel 2005 ) 
catalogs which include stars with proper motions 
larger than 0.18 mas yr“^ and 0.15 mas yr“^, 
respectively. We retained only the stars with a 
relative proper motion measurement error of less 
than 30%. The large and precise proper motions 
of these stars are useful to limit the number of false 
positives, i.e., common proper motion between two 
sources resulting from chance or measurement er¬ 
rors. To calculate the proper motion of the sources 
in the vicinity of these target stars , we c ompared 
the 2MASS PSC data ( Cutri et al. 20 0^ to data 


from the SDSS (Ac elman-McCa rthv et al. 2008! ) 


and/or the SIMP ( Artigau et al.l 20091) surveys. 
Of the roughly 60,000 stars included in the NLTT 
and LSPM catalogs, about 25,000 have accu¬ 
rate proper motion measurements and were ob¬ 
served by SIMP or SDSS; those constitute our 
initial sample. Target baseline’s between the 
SDSS/SIMP and 2MASS observations span from 
1 to 11 years, and the typical uncertainty of the 
calculated proper motions was 40 mas yr“^. 


2.2. Identifying Companions 


For each target star, we calculated the proper 
motion of all sources within a radius of 4' and 
selected any source with a proper motion within 
40 mas yr“^ of the target star as a candidate com¬ 
panion. Visual inspection was performed for all 
candidates to reject false detections of artifacts 
that are identified within catalogs as point sources. 
The available photometry from the 2MASS PSC 
catalog was then used to estimate a spectral type 
for each candidate companion (and target star 
when unknown). We also estimated a photomeric 
distance and its uncertainty based on the J-band 
magnit ude and the spectra l type-magnitude rela¬ 
tion of lHawlev et al.l ( 20021) when the information 
was not available in the literature. We required 
that the estimated distances of the target and can- 
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Table 1 

Photometric data 


Name 

2MASS Name 


a 


5^ 


^'b 

^/b 




K 

NLTT 251^ 

2MASS J00064746-0852350 

00^06^^ 

’^47'® .47 

-08 “ 

52' 

35.l" 

17.30 ± 0.01 

14.94 ± 0.01 

13.73 ± 0.01 

11.97 ± 0.02 

11.43 ± 0.02 

11.09 ± 0.02 


2MASS J00064916-0852457 

00^06^ 

’^49'® .17 

-08 “ 

52' 

45.7" 

20.91 ± 0.05 

18.26 ± 0.01 

16.45 ± 0.01 

14.14 ± 0.03 

13.55 ± 0.03 

13.13 ± 0.04 

NLTT 687 

2MASS J00140212-1814578 

00^14'^ 

'^02^.10 

-18 ° 

14' 

58.0" 




10.38 ± 0.02 

9.81 ± 0.02 

9.57 ± 0.02 


2MASS J00135882-1816462 

00^13^ 

'^58® .80 

-18 ° 

16' 

46.2" 




16.54 ± 0.14 

15.89 ± 0.18 

15.04 ± 0.13 

NLTT 2274*^ 

2MASS J00415543-I-1341162 

00^41^ 

’^55^ .43 

+ 13 ° 

4l' 

16.2" 

. . . g 

15.42 ± 0.01 

. . . g 

10.16 ± 0.03 

9.57 ± 0.03 

9.35 ± 0.02 


2MASS J00415453-I-1341351 

00^41'^ 

■^54^ .53 

+ 13 ° 

4l' 

35.l" 

21.15 ± 0.07 

18.71 ± 0.01 

16.83 ± 0.01 

14.45 ± 0.03 

13.67 ± 0.04 

13.24 ± 0.03 

BD-06 813 

2MASS J04050271-0600261 

04^^05^ 

'^02'® .71 

-06 “ 

oo' 

26.l" 

. . . g 

... g 

. . . g 

8.35 ± 0.02 

7.93 ± 0.03 

7.87 ± 0.03 


2MASS J04050209-0600409 

04^05^^ 

’^02'®.10 

-06 “ 

oo' 

40.9" 

21.22 ± 0.06 

18.54 ± 0.01 

17.08 ± 0.01 

15.11 ± 0.05 

14.30 ± 0.05 

13.88 ± 0.05 

NLTT 20640® 

2MASS J08583671 + 2711068 

08^58^^ 

’^36'® .72 

+ 27 ° 

ll' 

06.8" 

15.92 ± 0.01 

... g 

13.50 ± 0.01 

11.99 ± 0.02 

11.40 ± 0.02 

11.10 ± 0.01 


2MASS J08583693-f2710518 

08^58^ 

'^36® .94 

+ 27 ° 

lO' 

51.8" 

22.05 ± 0.11 

19.45 ± 0.02 

17.64 ± 0.02 

15.05 ± 0.05 

14.23 ± 0.05 

13.66 ± 0.05 

LSPM J1021+3704 

2MASS J10215240-H3704289 

10^21'^ 

'^52® .40 

+ 37 ° 

04' 

29.0" 

16.07 ± 0.01 

... g 

13.87 ± 0.01 

12.43 ± 0.02 

11.95 ± 0.02 

11.76 ± 0.02 


2MASS J10215386-(-3704166 

10^21'^ 

'^53® .87 

+ 37° 

04' 

16.7" 

24.67 ± 0.48 

21.85 ± 0.14 

19.91 ± 0.08 

17.18 ± 0.03^ 

16.25 ± 0.22 

15.66 ± 0.23 


2MASS J10432398-1706024 

10^43'^ 

'^23® .99 

-17 ° 

06' 

02.4" 




11.48 ± 0.02 

10.89 ± 0.02 

10.60 ± 0.03 


2MASS J10432513-1706065 

10^43'^ 

'^25® .13 

-17 ° 

06' 

06.6" 




15.67 ± 0.06 

15.07 ± 0.06 

14.61 ± 0.11 

NLTT 26746 

2MASS J11150134-t-1606447 

11^15'^ 

'^01® .35 

+ 16 ° 

06' 

44.7" 

. . . g 

16.61 ± 0.01 

... g 

11.14 ± 0.02 

10.56 ± 0.02 

10.32 ± 0.02 


2MASS 311150150+1607026 

11^15^ 

'^01® .50 

+ 16 ° 

;07' 

02.7" 

22.94 ± 0.37 

20.71 ± 0.07 

18.59 ± 0.05 

16.40 ± 0.12 

15.22 ± 0.09 

14.56 ± 0.11 

NLTT 29392 

2MASS J12024963+4204475 

12^02'^ 

'^49® .63 

+42 ° 

04' 

47.5" 

18.21 ± 0.01 

... g 

13.87 ± 0.01 

12.43 ± 0.02 

11.95 ± 0.02 

11.67 ± 0.02 


2MASS J12025009+4204531 

12^02'^ 

■^50® .09 

+ 42 ° 

04' 

53.2" 

21.45 ± 0.07 

18.87 ± 0.01 

16.96 ± 0.01 

14.50 ± 0.03 

13.69 ± 0.03 

13.26 ± 0.03 

LSPM J12594-1001 

2MASS J12594217+1001407 

12^59'^ 

'^42®.17 

+ 10 ° 

Ol' 

40.7" 

17.76 ± 0.01 

15.63 ± 0.00 

14.44 ± 0.00 

12.60 ± 0.02 

12.03 ± 0.03 

11.63 ± 0.02 


2MASS J12594167+1001380 

12^59^^ 

"^41® .68 

+ 10 ° 

Ol' 

38.06" 

24.45 ± 0.76 

21.58 ± 0.13 

19.60 ± 0.10 

16.72 ± 0.19 

15.53 ± 0.14 

14.89 ± 0.15 

NLTT 36369 

2MASS J14081918+3708294 

14^08'^ 

'^19®. 19 

+ 37° 

08' 

29.4" 

15.60 ± 0.01 

14.53 ± 0.01 

13.48 ± 0.01 

12.06 ± 0.02 

11.45 ± 0.02 

11.22 ± 0.02 


2MASS J14081969 + 3708255 

14^08’^ 

"^19®.70 

+ 37° 

08' 

25.6" 

21.98 ± 0.10 

19.09 ± 0.02 

17.47 ± 0.02 

15.42 ± 0.05 

14.78 ± 0.06 

14.35 ± 0.08 

LSPM J1441 + 1856 

2MASS J14412209 + 1856451 

14^41'^ 

■^22® .09 

+ 18 ° 

56' 

45.l" 

17.43 ± 0.04 

... g 

. . . g 

11.44 ± 0.02 

10.82 ± 0.03 

10.55 ± 0.02 


2MASS J14412565 + 1856484 

14^41'^ 

■^25® .65 

+ 18 ° 

56' 

48.5" 

23.38 ± 0.27 

21.34 ± 0.07 

19.41 ± 0.06 

16.94 ± 0.18 

16.17 ± 0.17 

15.40 ± 0.13 

NLTT 41701 

2MASS J15590815 + 3735477 

15^59'^ 

■^08® .15 

+ 37° 

35' 

47.7" 

. . . g 

... g 

. . . g 

9.12 ± 0.02 

8.67 ± 0.02 

8.56 ± 0.02 


2MASS J15590740 + 3735275 

15^59^ 

’^07® .40 

+ 37° 

35' 

27.5" 

22.89 ± 0.27 

19.66 ± 0.03 

18.07 ± 0.03 

15.75 ± 0.07 

14.92 ± 0.07 

14.53 ± 0.08 

HD 234344 

2MASS J16460765 + 5020405 

16'^46'^ 

"^07® .64 

+ 50 ° 

20' 

41.l" 

. . . g 

... g 

. . . g 

7.62 ± 0.02 

7.06 ± 0.02 

6.97 ± 0.03 


2MASS J16461148 + 5019456 

16^46^ 

"^ll® .48 

+ 50 ° 

19' 

45.7" 

20.47 ± 0.03 

17.43 ± 0.01 

16.66 ± 0.01 

13.61 ± 0.02 

12.98 ± 0.03 

12.59 ± 0.02 

HD 217246 

2MASS J22591790 + 0806484 

22^^59^ 

"^17® .90 

+ 08 ° 

06' 

48.4" 

. . . g 

... g 

. . . g 

8.56 ± 0.02 

8.26 ± 0.03 

8.19 ± 0.02 


2MASS J22591631 + 0806556 

22^59^^ 

’^16®.31 

+ 08 ° 

06' 

55.6" 

21.83 ± 0.09 

18.84 ± 0.01 

17.61 ± 0.02 

15.76 ± 0.06 

15.14 ± 0.08 

14.74 ± 0.10 

NLTT 56936 

2MASS J23274840+0451241 

23^27'^ 

'^48® .40 

+ 04 ° 

5l' 

24.2" 

. . . g 

... g 

. . . g 

8.16 ± 0.03 

7.58 ± 0.05 

7.41 ± 0.02 


2MASS J23274947+0450583 

23^27^ 

'^49® .48 

+ 04 ° 

50' 

58.4" 

21.54 ± 0.125 

18.92 ± 0.02 

17.17 ± 0.02 

15.10 ± 0.03 

14.36 ± 0.04 

13.97 ± 0.05 

TYC 1725-344-1 

2MASS J23553113+1755239 

23'^55m31'® .14 

+ 17° 

55' 

23.9" 

. . . g 

... g 

. . . g 

8.80 ± 0.03 

8.38 ± 0.02 

8.28 ± 0.02 


2MASS J2355345 + 175404 

23^55’' 

'^34® .50 

+ 17 ° 

54' 

04.0" 

23.24 ± 0.33 

20.43 ± 0.04 

18.73 ± 0.03 

16.39 ± 0.12 

15.42 ± 0.15 

15.02 ± 0.15 

NLTT 182 

2MASS J00054143+0626256 

00^05^ 

'^41® .43 

+ 06 ° 

26' 

25.6" 

17.08 ± 0.01 

15.42 ± 0.04 

14.51 ± 0.01 

13.02 ± 0.03 

12.51 ± 0.03 

12.19 ± 0.03 


2MASS J00054171 + 0626300 

00^05’^ 

"^41® .71 

+ 06 ° 

26' 

30.l" 

22.38 ± 0.16 

19.86 ± 0.03 

18.00 ± 0.02 

15.81 ± 0.10 

15.13 ± 0.11 

14.62 ± 0.11 

HD 2292 

2MASS J00265848 + 1705088 

00^26'^ 

'^58® .49 

+ 17° 

05' 

08.8" 

8.88 ± 0.01 

... g 

. . . g 

7.90 ± 0.027 

7.55 ± 0.04 

7.50 ± 0.03 


2MASS J00265989 + 1704463 

00^26^ 

"^59® .90 

+ 17 ° 

04' 

46.4" 

21.82 ± 0.13 

19.16 ± 0.02 

17.39 ± 0.02 

15.40 ± 0.05 

14.67 ± 0.05 

14.39 ± 0.07 

NLTT 4558 

2MASS J01221872 + 0330470 

01^22^ 

'^18®. 72 

+ 03 ° 

30' 

;47.l" 

. . . g 

... g 

. . . g 

7.42 ± 0.02 

7.14 ± 0.03 

7.07 ± 0.02 


2MASS J01221697+0331235 

01^22^^ 

’^16®.98 

+ 03 ° 

3l' 

23.6" 

22.05 ± 0.13 

19.95 ± 0.03 

17.95 ± 0.03 

15.47 ± 0.05 

14.65 ± 0.06 

14.40 ± 0.08 

NLTT 30510 

2MASS J12221994+3643539 

12^22'^ 

'^19®.94 

+ 36 ° 

43' 

54.0" 

. . . g 

... g 

. . . g 

10.51 ± 0.02 

9.98 ± 0.02 

9.73 ± 0.02 


2MASS J12221837+3643485 

12^22'^ 

■^18®.37 

+ 36 ° 

43' 

48.5" 

22.61 ± 0.17 

20.21 ± 0.03 

18.15 ± 0.03 

15.97 ± 0.08 

15.27 ± 0.09 

14.85 ± 0.09 

LSPM J1236-I-3000 

2MASS J12363558 + 3000315 

12^36'^ 

"^35® .59 

+ 30 ° 

oo' 

31.6" 

19.63 ± 0.02 

17.75 ± 0.01 

16.81 ± 0.01 

15.28 ± 0.05 

14.60 ± 0.06 

14.30 ± 0.08 


2MASS J12363647+3000315 

12^36^ 

’^36® .48 

+ 30 ° 

oo' 

31.5" 

23.04 ± 0.31 

20.63 ± 0.05 

18.86 ± 0.05 

16.82 ± 0.20 

16.06 ± 0.22 

15.48 ± 0.20 

NLTT 33793^ 

2MASS J13205010 + 0955582 

13^20^ 

’^50® .11 

+ 09 ° 

55' 

58.3" 

. . . g 

... g 

. . . g 

7.89 ± 0.03 

7.32 ± 0.04 

7.22 ± 0.03 


2MASS J13204159 + 0957506 

13^20^ 

"^41® .59 

+ 09 ° 

57' 

;50.6" 

20.40 ± 0.03 

17.47 ± 0.01 

15.80 ± 0.01 

13.73 ± 0.03 

13.08 ± 0.03 

12.61 ± 0.03 


^Data are from the 2MASS catalog ICJutri et al,i r2UUa) 

^Data are from the SDSS catalog <Ahn et aljrm^ 

^This system has been identifified to be a hierarchical triple system with an M6 primary and an M94-T5 secondary by [Burgasser_^t_^l ]Gni3) 

^This system has been identifified as a binary by [Fahert^_et_a^ J2Uiy) 

®This system has been identifified as a binary by [Zhang_et_a^ J2U1^ 

^This magnitude has been extrated from a CPAPIR MKO J-band image and has been converted to the 2MASS magnitude system. This object was poorly detected in the 2MASS calatog as it is 
faint. 
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didate companion agree within the estimated un¬ 
certainties. Following this procedure, we identified 
29 pairs of potential binary systems. We obtained 
imaging and spectroscopic follow-up observations 
of most of these candidates to improve their proper 
motion measurements and their spectral type de¬ 
terminations. We then applied a more detailed 
analysis, described below, to make a better assess¬ 
ment of their nature. Following this analysis, we 
found that 6 candidates systems are not comoving 
while 23 systems were found to be likely bound. 
These are the systems to be discussed in this paper 
(see Table [T]). 


Five of the 23 candidate pairs were found by 
other teams and reported as binary systems after 
we performed our initial search. The primary 
of these systems are NLTT 2274 (M4-I-M9.5), 
and NLTT 33793 (K4-HVt7.5), bot h found by 
Fahertv et al.l ( 2010( 1. NLT T 20640 (IZh ang et ;il1 


201(1 M4-^L01. NLTT 251 (iBurgasser et all II 2 OI 2 L 


M7-I-M9-I-T5), and HD234344 (IMason et al.ll2001a , 


K7,M7). In addition, shortly prior to the sub- 

mission of t his _ manuscript, a new study by 

Deacon et al. ( 2014[ l reported the discovery of 4 


other of our systems : NLTT 4558 (G5-I-L1-I-T3), 
NLTT56936 (K5-kM8), NLTT26746 (M4-HL4) and 
NLTT30510 (M2-HM9.5). 


3. NEAR-INFRARED IMAGING 


3.1. Observations 


We obtained J-band observations of 14 can¬ 
didate companions between 2012 November and 
2013 May using the Observatoire du Mont Megantic 
(OMM~) wide field n ear-infrared camera, CPAPIR 
( Artigau et al. 2004 1. installed on the OMM 1.6 m 
telescope, and t he OMM nea r-infrared spectro- 
imager SIMON ( Albert 2006ll installed on the 
CTIO 1.5 m telescope. The field of view (FOV) 
of CPAPIR at the OMM is 30'x30' with a pixel 
scale of 0.89" pixel“^, while the FOV of SIMON 
at the Cerro Tololo Inter-American Observatory 
(CTIO) 1.5 m is 8.4'x8.4' with a pixel scale of 
0.49" pixel“^. Typically, the observations con¬ 
sisted of 15 individual exposures of 3 to 30 s, with 
a dither step of 2' between them. These obser¬ 
vations provide a total time baseline of up to 15 
years with respect to 2MASS for proper motion 
measurements. 


3.2. Data Reduction 

The reduction of the imaging data followed 
standard procedures. The images of SIMON and 
CPAPIR were sky subtracted using a median sky 
image constructed from the entire data set taken 
over the night in the J filter. The images were flat- 
fielded using a fiat image derived from on and ojf 
dome images taken at the beginning or the end of 
the nights. Individual images astrometry was per¬ 
formed by cross-correlating point sources in the 
field with the 2MASS catalog. All images were 
then median-combined into a single science frame 
that was used for astrometric measurements. 

3.3. Results 

We measured the positions of the sources in our 
CPAPIR and SIMON follow-up imaging by finding 
all of the stars in the image and then fitting a 2D 
Gaussian over each one. These positions have been 
compared to those of all the objects in the 2MASS 
catalog that are up to 10' away from the target and 
the corresponding sources were associated. We 
determined the astrometric errors by computing 
the standard deviation of the difference between 
the positions in the CPAPIR/2MASS images and 
the corresponding positions in the 2MASS cata¬ 
log propagated to the CPAPIR/SIMON epoch ac¬ 
cording to the proper motion listed in the NO¬ 
MAD catalog (IZacharias et al.l l2004ll : this was 
done for reference stars of brightness comparable 
to our targets. We determined the source positions 
in the SIMP data, also obtained with CPAPIR, 
in the same manner; those measurements are also 
given in Table [5] as they have not yet been pub¬ 
lished elsewhere. 

For each candidate companion, we combined 
the position measurements from our follow-up 
imaging, 2MASS, SDSS, SIMP, and WISE to com¬ 
pute their proper motion using a weighted linear 
regression between the positions of the object and 
the measurement epochs. We adopted the uncer¬ 
tainty of the slope parameter as our uncertainty 
on the proper motion. The proper motions for 
all of our candidates, corrected for the parallax 
estimated from their photometric distance (see 
below), are compiled in Table [5] As an example, 
our calculated proper motions for all stars within 
10' of NLTT 26746 are plotted on Figure [1] ; the 
common proper motion of our candidate compan- 
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Table 2 

Astrometric measurements based on our observations 


Name 

RA 

Uncertainty ii 

ti RA 

DEC 

Uncertainty in DEC 

Epoch 

Source 


(deg) 

(mas) 


(deg) 

(mas) 

(Reduced JD) 


2M1043-1706A 

160.84974 

340 


-17.10118 

400 

54168.75 

SIMP 


160.84956 

560 


-17.10128 

670 

56345.70 

SIMON 


160.84953 

130 


-17.10125 

130 

56259.97 

CPAPIR 

2M1043-1706B 

160.85452 

340 


-17.10234 

400 

54168.75 

SIMP 


160.85435 

220 


-17.10245 

110 

56345.70 

SIMON 


160.85433 

250 


-17.10240 

130 

56259.97 

CPAPIR 

2M1646 + 5019 

251.54744 

190 


50.33065 

190 

55232.88 

SIMP 


251.54735 

400 


50.33103 

180 

56330.86 

CPAPIR 

LSPM12594-1001 

194.92540 

860 


10.02808 

720 

55232.91 

SIMP 


194.92522 

680 


10.02804 

270 

56330.82 

CPAPIR 

2M1259 + 1001 

194.92335 

860 


10.02741 

720 

55232.91 

SIMP 


194.92339 

320 


10.02728 

220 

53900.55 

CPAPIR 


194.92327 

380 


10.02737 

270 

56330.82 

CPAPIR 

LSPMJ1441 + 1856 

220.34193 

380 


18.94546 

390 

55232.92 

SIMP 


220.34189 

160 


18.94529 

130 

56330.89 

CPAPIR 

2M1441 + 1856 

220.35673 

380 


18.94627 

390 

55232.92 

SIMP 


220.37645 

260 


19.09413 

290 

56330.89 

CPAPIR 

NLTT182 

1.42300 

150 


6.44028 

170 

53980.73 

SIMP 

2M0005+0626 

1.42296 

150 


6.44028 

170 

53980.73 

SIMP 


1.43611 

360 


6.51388 

410 

56619.50 

SIMON 

NLTT20640 

134.65331 

310 


27.18477 

350 

55144.94 

SIMP 


134.65344 

120 


27.18451 

170 

56329.48 

CPAPIR 

2M0858 + 2710 

134.65424 

310 


27.18052 

350 

55144.94 

SIMP 


134.67524 

300 


27.31993 

340 

56329.48 

CPAPIR 

NLTT33793 

200.26834 

190 


10.02453 

180 

53838.67 

SIMP 


200.26833 

360 


10.02452 

280 

56330.83 

CPAPIR 

2M1320+0957 

200.17290 

190 


9.96377 

180 

53838.67 

SIMP 


200.17247 

110 


9.96350 

90 

56330.83 

CPAPIR 

NLTT36369 

212.07912 

360 


37.14177 

250 

55232.91 

SIMP 


212.08075 

330 


37.14089 

870 

56330.80 

CPAPIR 

2M1408 + 3708 

212.07911 

360 


37.14177 

250 

55232.91 

SIMP 


212.08084 

380 


37.14090 

390 

56330.80 

CPAPIR 

NLTT41701 

239.78411 

470 


37.59615 

400 

55232.89 

SIMP 


239.78417 

150 


37.59600 

160 

56430.58 

CPAPIR 

2M1559+3735 

239.78101 

470 


37.59051 

400 

55232.89 

SIMP 


239.78105 

120 


37.59038 

100 

56430.58 

CPAPIR 

2M0122+0331 

20.57078 

210 


3.52294 

150 

54036.68 

SIMP 

NLTT56936 

352.02213 

130 


4.97526 

170 

53986.70 

SIMP 


351.95329 

290 


4.85733 

250 

56572.54 

SIMON 

2M2327+0450 

351.95691 

130 


4.84984 

170 

53986.70 

SIMP 


351.95778 

70 


4.85024 

290 

56572.54 

SIMON 

2M0026+1704 

6.74936 

220 


17.07940 

280 

53994.73 

SIMP 

BD-06813 

61.26145 

470 


-6.00778 

470 

56349.54 

SIMON 


61.26139 

140 


-6.00779 

190 

56258.78 

CPAPIR 

2M0405-0600 

61.25892 

180 


-6.01190 

180 

56349.54 

SIMON 


61.25890 

110 


-6.01187 

110 

56258.78 

CPAPIR 

NLTT2274 

10.48038 

360 


13.68728 

220 

56619.54 

SIMON 

2M0041 + 1341 

10.47664 

140 


13.69254 

70 

56619.54 

SIMON 

2M2355+1754 

358.89345 

360 


17.90071 

140 

56619.54 

SIMON 

NLTT687 

3.50881 

1080 


-18.24962 

2160 

56547.70 

SIMON 

2M0013-1816 

3.49509 

1080 


-18.27954 

2160 

56547.70 

SIMON 

2M2259+0806 

344.81841 

140 


8.11571 

140 

56570.54 

SIMON 

2M0006-0852 

1.70457 

250 


-8.88072 

220 

56547.75 

SIMON 

LSPMJ1021 + 3704 

155.46781 

310 


37.07414 

70 

56329.67 

CPAPIR 

2M1021 + 3704 

155.47389 

490 


37.07081 

190 

56329.67 

CPAPIR 

NLTT26746 

168.75463 

130 


16.11187 

90 

56329.69 

CPAPIR 

2M1115 + 1607 

168.75524 

220 


16.11687 

150 

56329.69 

CPAPIR 

NLTT30510 

185.58421 

180 


36.73147 

110 

56329.73 

CPAPIR 

2M1222 + 3643 

185.57769 

290 


36.72998 

300 

56329.73 

CPAPIR 

LSPMJ1236+3000 

189.14888 

120 


30.00824 

110 

56329.71 

CPAPIR 

2M1236 + 3000 

189.15272 

410 


30.00832 

210 

56329.71 

CPAPIR 

2M1202+4204 

180.70753 

580 


42.08037 

640 

56330.73 

CPAPIR 
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ion with the primary is clearly seen. 



Fig. 1.— Proper motions for all objects within 
10' of NLTT 26746 based on our follow-up imag¬ 
ing (2MASS, SDSS, CPAPIR, and WISE data). 
The proper motion of the candidate companion 
2M1115-I-1607 is shown with a red diamond, while 
the proper motion of the primary star is shown 
with a red square. The two components of the 
candidate binary have proper motions consistent 
within the estimated uncertainties. 


4. SPECTROSCOPY 
4.1. Optical spectroscopy 

We obtained optical spectroscopy follow-up 
observations of all candidates using 
Multi-Object Spectrographs (CMOS) 

I 2 OO 4 II at the Gemini South and North telescopes 
during semesters 2008A and 2008B (programs 
GN-2008A-Q-80, GN-2008A-Q-20, GS-2008A-Q- 
12, GS-2008B-Q-34, and GN-2008B-Q-107). We 
used a 1" wide slit, the R400 grating with OG515 
blocking filter, and 2-pixel binning in both the 
spatial and the spectral direction. For each target, 
we obtained exposures with three central wave¬ 
length settings, 790 nm, 800 nm and 810 nm, to 
cover the small gaps between the three GMOS 
detectors. The exposure times range from 10s to 
600s (see Table [3|). We did not move the target 
in the slit during the observations. The result¬ 
ing spectra have a resolving power of R ^^850 
and cover the wavelength range from 650 nm to 
1000 nm. The typical seeing during the observa¬ 
tions was 0.75"-1.0". Standard calibrations were 
obtained for each science data set. We used the 


the Gemini 
(Hook et al. 


calibration data of the G subdwarf LTT9239, the 
DA white dwarf EG 274 and the BO star HZ44 for 
telluric and instrumental transmission calibration. 
Our main goal was to characterize the candidate 
companions, but given the short additional time 
needed, we also acquired a spectrum of the pri¬ 
mary star during a visit to a given system. When 
the primary was not too bright to saturate the 
detector in an exposure of a few hundred sec¬ 
onds, both the primary and the companion were 
put in the slit simultaneously, otherwise a shorter 
observation of the primary was obtained after ob¬ 
servation of the companion. Due to various issues 
during the execution of the observations, we were 
unable to procure data for all of the systems. In 
all, we obtained a spectrum for 14 candidate com¬ 
panions out for 23 and 10 candidate primaries out 
23. A log of the observations is given in Tabled 
The GMOS spectroscopy data were reduced us¬ 
ing a custom Interactive Data Language [IDL) 
routine. First, we corrected the individual expo¬ 
sures for the bias and flat field. We then extracted 
the spectral trace using a Moffat function extrac¬ 
tion profile, allowing for a linear sky solution. We 
subsequently constructed a quadratic wavelength 
solution using calibration arc lamps, and corrected 
the spectra for wavelength-dependent instrumen¬ 
tal transmission using a standard spectrum. The 
reference spe c trosco pic s pectra were taken from 
Hamuv et al. (1994) and Massev et al. ( 1988h . In 


the last step, we combined the individual spectra 
extracted from each of the 3 detector section, and 
finally combined all exposures together for a given 
object. The measurement uncertainties were esti¬ 
mated by extracting an empty sky region near the 
spectral trace, using the same method. All these 
spectra are normalized at 0.75 /iin. Figure [2] shows 
the GMOS spectra for 6 candidate primaries for 
which the spectral types were not published in the 
literature before this study. GMOS spectra at R 
^ 850 of ten candidates are shown in Figure[3]and 
the remaining four are shown in Figured with R 
400. The latter figure present the smoothed spec¬ 
tra our four faintest targets observed with GMOS. 
In all cases, a residual signal caused by an fringing 
effect from the detector is present. However, this 
effect affects the quality of our data significantly 
only for our faintest targets but even with this 
fringing, we are still able to identify the spectral 
type of the objects. 
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Fig. 2.— CMOS spectra of our 6 candidate primaries for which a spectral type was not available in the 
literature, displayed with a resolution of R ^ 850. T he spectra are norma lized at 0.75 /iin. As a comparison, 
we display in red the SDSS template spectra from Hawlev et al.l ( 2002h for the appropriate spectral type 
convolved to the same resolution. The feature in our spectra at 0.76//m is telluric absorption that we did 
not correct. 
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Fig. 3.— CMOS spectra of 10 of our candidate companions, displayed with a resolution of R ^ 850. The 
spectra are n o rmali zed at 0.75 fim. As a comparison, we display in red the SDSS template spectra from 
Hawley et al. ( 2002h for the appropriate spectral type convolved to the same resolution. The feature in our 


spectra at 0.76/rm is telluric absorption that we did not correct. Different observation setups account for the 
different wavelength range. 





















X (|am) 


Fig. 4.— CMOS spectra of 4 of our candidate companions, displayed with a degraded resolution of R ^ 
400 as these candidates were fainter and a lower signal-to-noise ratio was reached at full resolution. There is 
fringing in the data leading to a significant residual signal for these faint targets. T he spectra are normal ized 
at 0.75 ^m. As a comparison, we display in red the SDSS template spectra from iHawlev et al.l (|2002ll for 
the appropriate spectral type. 
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Fig. 

1700. 

IRTF 



5.— GNIRS spectra of 3 of our candidate companions, displayed with a spectral resolution of R ~ 
The spectra are normalized at 1.1 um.As a comparison , we display in red the spectra from SPEX 
Spectral Library ( Cushing et al. 2005; Ravner et al. 2009ll for the appropriate spectral type. 
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Table 3 

CMOS AND GNIRS Spectroscopic Observations 


System 

Component 

Obs. Date 

Instrument 

Number of exposures 

Exp. time 

Standard 



YYYY-MM-DD 



(a) 


NLTT251 

Both 

2008-06-14 

GMOS-S 

6 

600. 

EG 274 

NLTT 687 

Primary 

2008-06-15 

GMOS-S 

6 

60.4 

EG 274 


Companion 

2008-06-15 

GMOS-S 

6 

600. 

EG 274 

NLTT 2274 

Primary 

2008-08-05 

GMOS-S 

6 

60.4 

LTT 9239 


Companion 

2008-08-05 

GMOS-S 

3 

900. 

LTT 9239 

BD-06 813 

Primary 

2008-01-28 

GMOS-S 

6 

60.4 

EG 274 


Companion 

2008-01-28 

GMOS-S 

6 

600. 

EG 274 

NLTT 20640 

Both 

2008-02-11 

GMOS-N 

6 

599. 

HZ44 

LSPM J1021-I-3704 

Both 

2008-02-27 

GMOS-N 

2 

599. 

HZ44 

2M1043-1706 

Primary 

2008-01-28 

GMOS-S 

6 

600. 

HZ44 

NLTT 26746 

Primary 

2009-01-01 

GMOS-S 

3 

900. 

LTT 9239 

NLTT 29392 

Both 

2008-02-11 

GMOS-N 

6 

599. 

HZ44 

NLTT 36369 

Both 

2008-04-14 

GMOS-N 

9 

599. 

HZ44 

LSPM J1441-f-1856 

Companion 

2007-08-06 

GMOS-S 

10 

600. 

LTT 9239 

NLTT 41701 

Primary 

2008-03-15 

GMOS-N 

6 

9.99 

HZ44 


Companion 

2008-03-15 

GMOS-N 

9 

599. 

HZ44 

HD 234344 

Primary 

2008-03-02 

GMOS-N 

6 

29.9 

HZ44 


Companion 

2008-03-02 

GMOS-N 

9 

419. 

HZ44 

HD 217246 

Companion 

2007-06-24 

GMOS-S 

6 

600. 

LTT9239 

NLTT 56936 

Companion 

2008-12-27 

GMOS-S 

3 

900. 

LTT9239 

TYC 1725-344-1 

Companion 

2007-06-21 

GMOS-S 

6 

600. 

LTT9239 

NLTT 182 

Primary 

2008-06-30 

GMOS-S 

6 

600. 

EG 274 

2M1043-1706 

Companion 

2007-04-01 

GNIRS 

8 

120. 

HD 137873 

2M1115-f 1607 

Companion 

2007-01-08 

GNIRS 

11 

60.0 

HD 137873 

2M1259-f 1001 

Companion 

2007-01-08 

GNIRS 

15 

60.0 

HD 137873 


4.2. Near-Infrared Spectroscopy 

Three of our candidate companions, 2M1043- 
1706B, 2M1115+1607, and 2M1259+1001, were 
identified on their own as brown dwarf candi¬ 
dates in a separate search for high proper motion 


brown dwarfs in the SIMP data (jArtigau et al 


2 OO 9 I Robert et ah, in preparation), but have not 


yet been published. As other candidates identi¬ 
fied in that search for brown dwarfs, they were 
followed-up using near-infrared spectroscopy with 
GNIRS at the Gemini South telescope in semester 
2007A (GS-2007A-Q-28). We used a 0.3" wide 
slit, the 32 1 mm“^ grating, the short camera, and 
the cross-dispersed mode for a resolving power of 
R ^ 1700 and a spectral coverage from 0.9 fim to 
2.5 /rm. Observatory standard calibrations were 
obtained with each observation. For telluric and 
instrumental transmission calibration, the AO star 
HD 137873 was observed shortly before/after each 
observation. The observations were obtained with 
a typical ABBA dither pattern along the slit with 
individual exposures of 60 to 120s with 8 to 15 ex¬ 
posures. The details of the observations are shown 


then corrected for both spectral and spatial dis¬ 
tortions. To extract the spectra, a positive and 
a negative extraction box that matched the trace 
profile was used. An argon arc lamp spectrum 
was taken at the end of each sequence, which was 
used as a first wavelength calibration. The wave¬ 
length scale was then mor e precisely adjusted t o 
the atmospheric OH lines (IRousselot et al.ll2000ll . 
Finally, the spectra for each target extracted from 
image pairs were normalized and their median was 
taken to make the final spectra. The GNIRS spec¬ 
tra of the three companions are shown on Figure [5] 
along with a s pectra from the SPEX IRTF Spec¬ 
tral Library ( Cushing et al.l 120051 : iRavner et al.l 
l2009li^ corresponding to the same spectral type as 
determined below. 

4.3. Results 

We determined the spectral type of all ob¬ 
jects observed with CMOS from spectral in¬ 
dices, and then confirmed them by visual com¬ 
parison with templates. For each spectrum we 
computed the following spectral indices: PCS, 


The reduction of the spectra was made using 
custom (IDL) routines as described in Delorme et al. 

ana iviarxin et ai. iiyyyi). iiuo ana vu irom 
Cruz & Reid (20^. CaH. Ti I. and Ca H from 

(2008ll. First, successive image pairs were sub- 

Kirkpatrick et al. (1991). and VO 2 and Ti07 de- 


tracted to remove the sky emission and divided 
by a median combined spectral flat. They were 


^http://irtfweb.ifa.hawaii.edu/-spex/ 
IRTF_Sp6Ctral_Library 
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fined bv fLepine et alJ ( 2003 '). Each of these indices 
constitutes a good spectral type indicator in some 
restricted spectral type interval, but altogether 
they cover the range from K5 to L6 dwarfs. The 
index values-spectral type relations and regimes of 
validity were taken from their respective papers. 
The measured indices are compiled in Table IH 
along with the corresponding spectral types. The 
latter are averages of the spectral types obtained 
from all the indices that are valid for a given type. 
The few blanks in the table are caused by the 
shorter spectral coverage of some spectroscopic 
observations. We visually compared ou r spec¬ 
tra t o templates from the Pickles Atlas ( PicklesI 


19981) or to SPSS sp ectra of M and L dwarfs from 


Hawley et al.l (120021) to confirm spectral typing. 


The spectral types of all companions observed 
with CMOS range from M6.5 to LI. 

We determined the spectral type of the three 
companions observed with GNIRS using th e spec¬ 
tral indices defined by Allers &: Liu ( 20131) . which 
are based on the H 2 O, H 2 O-I, H 2 O- 2 , and H 2 OD 
indices and which are valid for M5V to LTV 
dwarfs. Then, we visually compared our spectra to 
the templates of M and L dwarfs from the SPEX 
IRTF Spectral Library. We determined a spectral 
type of M9.0±0.5 for 2M1043-1706B, L5 ± 1 for 
2M1115-H1607, and L4.5 ± 0.5 for 2M1259-hl001; 
these are indicated in Tabled 

There are 19 objects for which we do not have 
CMOS nor GNIRS observations. Among these are 
10 primaries and 2 companions for which a spec¬ 
tral type was known from the literature and we 
adopted those values. For the remaining seven, we 
estimated a spectral type from their 2MASS and 
SDSS col ors based on the color-spect ral type re¬ 
lations of Sheppard fc Cushind ( 20091) . The spec¬ 
tral types adopted for all objects are summarized 
in Tabled This table is divided in the two parts: 
the top half is for systems for which we have a 
spectrum for the companions while the lower half 
is for systems for which we did not have a spec¬ 
trum for the companions. 

For each object we estimated a photomet¬ 
ric distance using their 2MASS J-band magni- 
tude and the Afj- spectral type relations from 
Hawley et al. (I2OO2I) for la te-K to M5 dwarfs, and 


from Dupuv fc Liu ( 20121) for later type objects, 
M6 to L5. These distance estimates are given in 
Tabled 


5. PROBABILITY OF RANDOM ALIGN¬ 
MENT 


After determining the proper motion, the spec¬ 
tral types and the distance for all components, we 
need to confirm that they form a bound pair. The 
probability of random alignment is the probability 
that from chance, in a search like ours, we would 
find a physically unrelated companion star having 
the same proper motion and photometric distance 
as the primary, within our uncertainties, and that 
is separated by less than the observed separation 
of our candidate. This is given by the probability 
of finding two stars of the relevant spectral types 
within a given separation of each other and at the 
same distance, times the probability of finding two 
stars of the relevant spectral types with the same 
proper motion, given that they are at a similar 
location on the sky. 


For the former probability, as our search tar¬ 
geted a fixed set of potential primary stars, we 
need only calculate the probability of finding a 
companion star close to those primaries. All of our 
candidate companions have spectral types in the 
M6-L5 range, thus this is the relevant range for 
statistical calculations. From the spatial density 
of M6-M8 dwarfs of 2.2 x 10“^ pc~^ per /-band 
magni tude interval determined by Phan-Bao et al.l 
( 2003 ) and the spatial density of M8-L3.5 of 1.64 x 
10“^ pc~^ per J-band ma g nitud e interval deter¬ 
mined by Phan-Bao et al.l ( 20081) . we calculated 
the overall spatial density of late-M dwarfs to 
early-L dwarfs (M6 to L3.5). We used this density 
to calculate the number of such stars in a spherical 
shell of radius equal to the distance from the Sun 
to a candidate and thickness given by our distance 
estimate uncertainty. Then, from this number we 
calculated, the average number of stars in a sky 
projected disk of radius equal to the separation of 
our candidate binary, assuming an isotropic distri¬ 
bution of stars. We multiplied this number by 25 
000 , the number of primary stars targeted by our 
search, and finally calculated the corresponding 
probability of detecting at least one object within 
this total search area. 


For the latter probability, we used Monte Carlo 
calculations to determine the proper motion dis¬ 
tribution of stars of a given spectral type and at 
the sky position of the candidate system based on 
the observed Galactic space velocity (UVW) dis- 
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Table 4 

Spectral Indices 


Object 

PC3®^ 

PC4^ 

PC5'‘ 

Ti05^ 

vo*> 

CaH*^ 

Ti 

Ca II*^ 

V02‘^ 

Ti07‘^ 

Spectral Type ® 

NLTT251 

1.3695 

1.6816 

1.7827 

0.19484 

2.1168 

1.5603 

1.0704 

1.0764 

0.56026 

0.65387 

M6 ± 0.5 

2M0006-0852 

2.3435 

3.8938 

4.8902 

0.34909 

2.3607 

1.5009 

1.0326 

1.0677 

0.32328 

0.46659 

M9 ± 0.5 

NLTT687 

1.0306 

1.0248 

0.93255 

0.58126 

2.0021 

1.1914 

1.0622 

1.2380 

0.89734 

0.92856 

M3 ± 0.5 

2M0013-1816 

2.8484 

6.2637 

7.3279 

0.96672 

2.3180 

1.4236 

1.0617 

1.1923 

0.60007 

0.95069 

Ll ± 0.5 

NLTT2274 

1.1037 

1.1221 

1.0918 

0.40346 

1.9970 

1.2715 

1.0803 

1.1374 

0.77698 

0.85842 

M4 ± 0.5 

2M00414-1341 

2.7069 

4.6043 

5.7158 

0.68651 

2.5128 

1.3602 

1.0093 

1.0645 

0.36967 

0.62286 

M9± 0.5 

2M0405-0600 

1.6339 

2.3019 

2.4712 

0.26199 

2.2343 

1.3667 

1.1098 

1.1253 

0.42198 

0.57595 

M6 ± 1 

NLTT20640 

1.1217 

1.2361 


0.37900 

2.0244 

1.3055 

1.0679 

1.1134 

0.72927 

0.85573 

M4 ± 0.5 

2M0858-|'2710 

2.4959 

4.7033 


0.64835 

2.3871 

1.2406 

0.98143 

1.0637 

0.38629 

0.65786 

LO ± 1 

LSPMJ1021-I-3704 

1.0898 

1.1628 


0.38845 

1.9939 

1.2819 

1.0796 

1.1583 

0.75780 

0.88121 

M4 ± 0.5 

2M1021-I-3704 

2.7920 

5.3320 


1.0516 

0.46435 

0.92451 


1.0668 

0.60661 

0.66820 

L2 

2M1043-1706A 

1.0665 

1.1986 

1.1118 

0.41786 

1.9828 

1.2421 


1.1573 

0.77439 

0.87779 

M4 ± 0.5 

NLTT26746 

1.0896 

1.0664 

1.1495 

0.44003 

1.9733 

1.2758 

1.1094 

1.1651 

0.83645 

0.85031 

M4 ± 0.5 

NLTT29392 

1.3948 

2.0186 


0.21082 

2.1644 

1.4311 

1.0827 

1.0974 

0.49018 

0.65474 

M6 ± 0.5 

2M1202 + 4204 

2.7094 

5.3242 


0.66899 

2.4729 

1.3948 

0.96756 

1.0718 

0.39638 

0.70797 

LO ± 0.5 

NLTT36369 

1.0359 

1.0841 


0.42531 

2.0007 

1.2596 

1.0741 

1.1489 

0.80273 

0.89932 

M4 ± 0.5 

2M1408-I-3708 

1.6805 

2.5229 


0.22558 

2.3573 

1.4836 

1.0556 

1.1543 

0.37181 

0.57255 

M7 ± 0.5 

2M1441-I-1856 

2.7019 

4.2794 

2.7264 

0.88903 

2.2368 

1.5198 

0.75043 

1.2276 

0.50723 

0.76976 

Ll ± 1 

2M1559-I-3735 

1.8393 

2.9506 


0.27357 

2.4630 

1.3777 

1.0318 

1.1406 

0.33012 

0.55086 

M8 ± 0.5 

2M1646-I-5019 

2.0450 

3.0980 


0.17044 

2.1751 

1.8641 

1.0559 

1.1482 

0.39253 

0.54864 

M7 ± 0.5 

2M2259-I-0806 

1.5610 

1.8795 

1.6930 

0.14590 

2.1279 

1.5757 

1.0905 

1.2289 

0.50888 

0.61669 

M6 ± 0.5 

2M2327 + 0405 

2.1572 

3.2263 

3.8603 

0.27870 

2.4432 

1.4433 

1.0444 

1.0839 

0.33743 

0.49614 

M9 ± 0.5 

2M23554-1754 

2.2927 

3.7816 

2.9846 

0.42313 

2.4953 

1.3452 

1.0084 

1.1327 

0.40961 

0.76501 

M9 ± 0.5 

NLTT182 

1.1426 

1.2647 

1.1793 

0.34680 

2.0167 

1.3121 

1.0885 

1.1310 

0.69285 

0.79268 

M4.5 ± .5 


^The spectral indice is from [^lartir^_et_alj 
^The spectral indice is from [r^ruz_^_Jie^ 

^The spectral indice is from[Rirk£atrick_et_j^ J^991D 
“^The spectral indice is from [rie£ine_et_^^ J2UU^ 

®Tlie spectral types shown in the last column are the median of the spectral types computed using the spectral indices presented in the table. 


tributions of stars in the solar neighborhood. The 
UVW distributions were derived from the v elocity 
dispersions given by Mi|rami_&_JIecy ([WS^ for F, 
G and K stars , by Bochanski et ai r ( 2011 ) for M 
dwarfs, and by Schmidt et al.l ( 201011 for L dwarfs. 
We then determined the probability that a star of 
the spectral class of the primary and a star of the 
spectral class of the secondary, at the sky position 
of the system, would have proper motions larger 
than 0.1 mas yr“^ and consistent with each other 
within 2(7. 


We finally multiplied the above two probabil¬ 
ities to get the probability of random alignment; 
the results are given in Table [S] For all of our 
candidate systems, this probability is less than 
1.3 X 10“^, indicating that they most likely form 
a physical pairs. 


6. DISCUSSION 

We estimated the mass of the companions in 
our systems based on evolution models; this re¬ 
quires first estimating an age for the systems. The 
tangential velocity of our systems, which is an in¬ 
dicator of the population to which an object be¬ 
longs, has been compared to the Besangon Galac¬ 


tic model ( Robin et al. 2003ll . All of our systems 
have tangential velocities consistent with being 
part of the thin disk. According to the Galac¬ 
tic model simulations, most of our targets have 
a > 80% probability to be younger than 3 to 5 
Gyr. We also looked for sign of youth in our 
GNIRS spectra. We used t he spectral indices de¬ 
fined by Allers Liu ( 2013h . already used to com¬ 
pute the spectral types, which are also sensitive to 
surface gravity. According to these indices, none 
of the three companions have low surface grav¬ 
ity. We then looked for Ha emission, at 656.3 nm, 
which is indicative of the st ellar activity and mod- 
erately correlated with age (|West et al.ll2008l) . Vi¬ 
sual inspection of the spectra and calculations of 
the Ha equivalent widths indicate that 2M0405- 
0600 (M6.5), 2M1202-b4204 (LO) and 2M1043B 
(M9) have equi valent widths higher t han 1 A, and 
are thus active. Schmidt et al.l ( 2010[) had already 
identified that NLTT 29392B was an active LO 
dwarf. However, this does not mean that they are 
particularly young as M dwarf s can be act i ve du r- 
ing a few Gyr (see discussion in West et ^ ( 20081) 1 


^See http://www.vlinbiiiaries.org/ for an up-to-date cen¬ 
sus of the known very low mass binaries, maintained by N. 
Siegler, C. Gelino and A. Burgasser. 
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Table 5 

Derived parameters for the binary systems 


Name 

Spectral type®" 

cos S 
(mas/yr) 

^5 

(mas/yr) 

Separation 

(") 

Distance^ 

(pc) 

Prob ® 

P. sep. 
(AU) 

Pos. angle 
(deg) 

Mass 

(Mq) 

NLTT 25lP 

2M0006-0852 

M6.0 ± 0.5 

M9.0 ± 0.5 

-70 ± 9 

-59 ± 18 

-311 ± 6 

-318 ± 9 

27.4 ± 0.6 

25 + "^ 

37±5 

-11 

2.8 Xl0“3 

83o1?52 

113 ± 1 

0.102-0.133 

0.079-0.085 

NLTT 687 

2M0013-1816 

M3.0 ± 0.5 

Ll.O ± 0.5 

-44 ± 12 

-20 ± 41 

-173 ± 5 

-212 ± 37 

118.1 ± 0.6 

42-t-.^D 

^^-16 

sstlj 

3.9 X10“^ 


203.1 ± 0.3 

0.389-0.412 

0.072-0.078 

NLTT 2274^ 

2M0041 + 1341 

M4.0 ± 0.5 

M9.5 ± 0.5 

-165 ± 56 

-162 ± 58 

-155 ± 2 

-153 ± 13 

23.2 ± 0.8 

23+9* 

40+3 

1.6 X10“^ 


325 ± 2 

0.207-0.272 

0.076-0.083 

BD-06 813 

2M0405-0600 

K0° 

M6.5 ± 1 

31 ± 8 

42 ± 6 

-139 ± 5 

-117 ± 10 

17.4 ± 0.2 

68 ± 10^ 

®®-37 

7.6 X10“'^ 

1340+^“^ 

-1-^4U-340 

211.6 ± 0.5 

0.925-1.014 

0.096-0.114 

NLTT 20640^ 

2M0858+2710 

M4.0 ± 0.5 

LO ± 1 

103 ± 5 

104 ± 7 

-179 ± 5 

-182 ± 6 

15.6 ± 0.2 

54+^^ 

°^-20 

2.2 X10“® 

^®°~270 

168.9 ± 0.6 

0.207-0.272 

0.074-0.081 

LSPM J1021+3704 

M4.0 ± 0.5 ^ 

-132 ± 9 

-141 ± 4 

22.2 ± 0.4 


1.3 Xl0“3 

30001675 

125 ± 1 

0.207-0.272 

2M1021 + 3704 

LO ± 1 

-131 ± 8 

-124 ± 12 






0.071-0.076 

2M1043-1706A 

2M1043-1706B 

M4.0 ± 0.5 

M9.0 ± 0.5 

-97 ± 16 

-90 ± 14 

-140 ± 22 

-141 ± 22 

17.1 ±0.1 


7.1 X10“^ 

102014“ 

84 ± 0.5 

0.207-0.272 

0.079-0.085 

NLTT 26746“ 

2M11154-1607 

M4.0 ± 0.5 

L5 ± 1 * 

-247 ± 8 

-253 ± 28 

-141 ± 2 

-127 ± 20 

18.0 ± 0.2 


1.5 X10“^ 

6601410 

6.5 ± 0.7 

0.207-0.272 

0.056-0.073 

NLTT 29392 

2M1202+4204 

M6.0 ± 0.5 

LO.O ± 0.5 

-306 ± 11 

-292 ± 6 

-300 ± 17 

-276 ± 11 

7.3 ± 0.1 

33l?4 

1.1 X10“® 

3ioi;55 

42.0 ± 0.8 

0.102-0.133 

0.074-0.081 

LSPM J1259+1001 

M5 ± iS 

-142 ± 15 

22 ± 9 

7.65 ± 0.08 

42"^^^ 

^*^-21 

2.9 X10“^ 

345+iyp 

■^40-156 

250.9 ± 0.6 

0.121-0.167 

2M1259+1001 

L4.5 ± 0.5 

-146 ± 5 

29 ± 7 


47±5 
^ -18 




0.057- 0.074 

NLTT 36369 

2M1408+3708 

M4.0 ± 0.5 

M7.5 ± 0.5 

-247 ± 8 

-249 ± 1 

80 ± 1 

80 ± 10 

7.9 ± 0.1 

00-21 

84+^ 

®^-29 

9.3 Xl0“3 

cqn+ lOii 
“^^^-180 

122 ± 1 

0.207-0.272 

0.087-0.096 

LSPM J1441 + 1856 

M6P 

-56 ± 9 

-154 ± 9 

51.1 ± 0.2 

OO-IO 

8.7 X10“^ 

4110+^^^^ 

4ttU-ii66 

86.4 ± 0.2 

0.102-0.133 

2M1441 + 1856 

LI ± 1 

-59 ± 9 

-185 ± 12 


qq±10 

^^-34 




0.072-0.079 

NLTT 41701 

2M1559+3735 

K2^ 

M8.0 ± 0.5 

42 ± 12 

46 ± 9 

-137 ± 9 

-152 ± 14 

22.4 ± 0.5 

89"^® 

®^-29 

8.1 X10“^ 

17351000 

204 ± 1 

0.872-0.941 

0.084-0.090 

HD 234344® 

2M1646+5019 

K7^ 

M7.5 ± 0.5 

-129 ± 29 

-118 ± 4 

394 ± 11 

400 ± 7 

69.3 ± 0.3 

31.6 ± 0.9"^ 

33l?4 

2.4 X10“^ 

270011“ 

146.2 ± 0.3 

0.705-0.751 

0.087-0.095 

HD217246 

2M2259+0806 

G5^ 

M6.5 ± 0.5 

137 ± 34 

122 ± 6 

59 ± 19 

69 ± 5 

24.7 ± 0.2 

{6S1^^)^ 

121 + 20 

3.6 X10“^ 

2340+^^^° 

-so-iU — j^QOO 

286.0 ± 0.5 

0.821-0.930 

0.096-0.113 

NLTT 56936“ 

K2-f-K5“ 

443 ± 1 

171 ± 8 

30.3 ± 0.1 

64 ± 9*^ 

4.7 Xl0“3 

1 825"*”^*“ 
lS25_3go 

147.8 ± 0.2 

1.46±0.09 “ 

2M2327+04505 

M9 ± 0.5 

443 ± 6 

191 ± 12 






0.079-0.085 

TYC 1725-344-1 

2M2355-I-1754 

G5III^ 

M9.0 ± 0.5 

-65 ± 21 

-69 ± 16 

-78 ± 14 

-86 ± 3 

93.8 ± 0.4 

103113 

2.8 X10“'^ 

0700ll“OO 

149.0 ± 0.23 

1.91±0.11 ^ 

0.079-0.085 

NLTT 182 

2M0005-I-0626 

M4.5 ± 0.5 

LO ± iS 

213 ± 6 

205 ± 7 

-113 ± 7 

-126 ± 8 

6.1 ± 0.1 


1.4 X10“^ 

400il» 

43 ± 4 

0.157-0.181 

0.079-0.085 

LSPM J1236-I-3000 

M6 ± iS 

140 ± 5 

-124 ± 4 

12.2 ± 0.8 


1.9 Xl0“3 

158011“ 

89 ± 4 

0.102-0.133 

2M1236-I-3000 

M9.0 ± iS 

162 ± 27 

-102 ± 15 


126ll| 




0.079-0.085 

HD 2292 

2M00264-1704 

G5^ 

M9 ± iS 

-143 ± 15 

-149 ± 2 

-113 ± 19 

-128 ± 41 

30.5 ± 0.5 

(42l4“)‘ 

65 + ^ 
^®-20 

2.6 X10“^ 

1200l“O 

138 ± 1 

1.05-1.12 

0.108-0.118 

NLTT 30510“ 

2M12224-3643 

M3 ± 1^ 

LO® 

269 ± 12 

272 ± 18 

-39 ± 16 

-36 ± 8 

20.7 ± 0.5 

^°-23 

70 ± 10^ 

1.9 X10“^ 

163510“ 

254.2 ± 0.9 

0.377-0.431 

0.074-0.081 

NLTT 33793J 

K5j 

-210 ± 64 

-131 ± 12 

168.5 ± 0.3 

(38.i12'3)J 

5.4 X10“® 

5780_ii9o 

311.7 ±0.1 

0.724-0.782 

2M13204-0957 

M8j 

-222 ± 38 

-141 ± 10 


36 ± 3^ 



0.087-0.096 

NLTT 4558“ 

G5“^ 

22 ± 13 

-142 ± 15 

44.8 ± 0.8 

(58 ± 3^ 

3.5 Xl0“3 

222213“ 

324.2 ± 1 

1.05-1.12 

2M0122-I-0331 

L2 ± 1® 

45 ± 10 

-162 ± 5 


43 + 13 

^•^-40 




0.071-0.076 


^■Spectral type are extracted from our CMOS or GNIRS spectra unless noted otherwise,^ Distance computed using the spectral types and the relations from 
IDupuv & cn 120121) and [Hawle^_et_a^ {2003)i ^ Probability of random alignement,^ Masses are evaluated using the mass—Mj relation from the BT-Settl model of 
lAllardI 120141) , TZhang er~an i2009al'l ,^ This system has been identified as a M4-f-L0 binary by [Fahert^_et_a^ (2013). Their spectral types, proper motions, distance 
and sepration are consistent with the ones we found,® Spectral type calculated from i’-J iSheppard &c Gushing||20 ^ IPickles fe Depagnl 120111) , * lAmmons et alj 
12003) has estimated a distance that matches to one obtained from the spectral type extracted from the optical GMOS spectrum, ® This system has been identified as 


a M4-f-L0 binary by| 


distance from 


et al. 


20^n,^ Heliocentric distance from ^ndersoi^_&_^]ranci3 (2013)^^_Photometric distance from^hang_^t_a^ (2009a^ 


2013 


Heliocentric 

_195^) ,P This system has been identified to be a hierarchical 
Their proper motions, distance and sepration are consistent 


_ _ _ (^993) 1 ° The spectr al type is from|^^rne^ 

triple system with an M7.0ih0.5 primary and an M8.5i0.5 + T5ihl secondary by [Burgasser^^^^yj^^^^. _^_^__,_____ 

with the ones we found,® This system has been identified as a M4+L0 binary by|2nan^^'^^rrf2^^5)^ Their spectral types, proper motions, distance and sepration 
are consistent with the ones we found,^ [Kharchenk^ (200^,® This system has been identified as a K7+M7 binary by [Mason_et_^ ] ji^UOlal) . Their spectral types, 
proper motions, distance and sepration are consistent with the ones we found^^| (Zhang et al.l (2QQ9b|) found a spectral type of LI by colors and has identified it as a 
brown dwarf candidate.,^ This system has been identified by [Deacor^^t_j^ (20^^",^ The mass of NLTT 4558 has been calculated by [Porto_de_^dello_^t_alJ (201^) 
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Fig. 6.— Separation (in AU) versus mass ratio for 
our systems and for known binaries. The red cir¬ 
cles represent our sample. The LSPMJ1259-I-1001 
system and the NLTT26746 system, for which 
the companion is either a L4.5 or a L5 dwarf, 
are shown with the open green circle. Very low 
mass binaries from the Very Low Mass Binaries 
Archiv(j^re shown as a diamond w hile the trian¬ 
gles r epresent binary systems from jGomes_et_aI 


(2012). Dhital et al.1 

(2010. Fahertv et al. (20inh 

and references therein, Raghavan et al. (20101). 

Metchev & Hillenbrandl ( 

2004), Reid et al. 

(200ll). Wilson et al. 

(2001) 

Mason et al. 

(2001b[l. 

Reid & GizisI (1997 

. Fischer & Marcv 

(19921), 


and Close et ^1 1 199dl . 


on the relatively weak constraint that Ha provides 
for an individual object). 


In addition, we used the BANYAN II code 
( Gagne et al.l l20l3l to test if any of our targets 
belong to a nearby young moving group. We 
found out that while most of our target have a 
higher probability of belonging to the field, one of 
them (NLTT182) has a 98.8% probability of be¬ 
ing a Beta Pictoris member, with a corresponding 
false positive rate of 4.5%. This is quite inter¬ 
esting but by no means a proof of youth or true 
membership; the complete 3D kinematics would 
have to be confirmed, as well as youth from var¬ 
ious indicators. We only have an optical spec¬ 
trum of the primary star of this system, and we 
looke d for gravity-sensitive spectroscopic indices 
in it JiC rkpatrick et al. 2000 : Cruz fc Reid 2002 : 


Cruz et al.boOT : Reid et ahl 2008 1. such as the Na 

a index, the Na-b index, the K-a index, and the 


CrH-a index, to assess whether this object is young 
or not. However, these four indices fall close to the 
field scatter and the result of our analysis is mainly 
inconclusive. More observations will be necessary 
to firmly establish whether or not this system is 
young. For the purpose of this paper, due to the 
lack of solid evidence, we simply assume this sys¬ 
tem to be of the age of the field, as we do for all 
of our other systems. 


Based on the above analyses, we thus esti¬ 
mated that the ages of all of the systems are 
roughly 1-7 Gyr. We then evaluated the masses 
using the mass-M.r re lation from the BT-Settl 
model of lAllar^ ( 2014 1. where the Mj are ob¬ 
tained fro m our inferred spect ral t ypes using the 
relatio n of lHawlev et al.1 (1200211 and lDupuv fc Liu 
()2012( l. In the cases where the distance were 
not well defined, we decided to use the effective 
temperature instead of the J magnitude to find 
the masses. The effective temperature for each 
target has been extracted fr om its spectral type 
accord ing to the relations in IPecaut fc Mamaiek 
( 2013l l. We note that masses infered in this man¬ 
ner are overestimated by about 10% for objects 
with a G or early-K spectr al types as compared to 


Henry fc McCarthy! (1199311 . The estimated masses 


for ages between 1 to 7 Gyr are given in Table [5] 
The uncertainties come from the propagation of 
the error of spectral types on the inferred magni¬ 
tudes (or effective temperatures). 

Figure [5] shows the separation in AU of our sys¬ 
tems as a function of their mass ratio, as com¬ 
pared with other known binaries. Some of our new 
systems, having low mass ratios, reach a sparsely 
populated regions of the diagram. In particular, 
we found 14 systems with a mass ratio less than 
0.3, below the bulk of the previously known pop¬ 
ulation. 


As some of our system have a low binding 
energy, it is interesting to find out if the sys¬ 
tems are stable. Orbital evolution is possible for 
wide, weakly bound binaries and it might lead to 
their disruption as they travel through the Galaxy 
and encounter stars a nd giants molecular clouds. 
Weinberg et al. ( 1987l l worked out the calculations 


for the evolution and lifetimes of such wide binary 
systems in the solar neighbourhood; these can eas- 
ily be scaled to very low-mass stars as was done by 


Artigau et al. (20^. For all of our binaries, 


we 

find that the half-life is of the order of a Hubble 
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time or larger. The probability of survival is thus 
high for all systems. 


Burgasser et al.l 


( 2005h have shown that the 
binary fraction is higher for ultra cold dwarfs 
that are in a wide binary system where the 
primary is a stellar o bj ect. Furthermore, 
IWhitworth fc StamatellosI ( 2006ll discussed the 
possibility that H 2 dissociation might trigger a 
secondary fragmentation of a companion low-mass 
protostar if the latter is at the cooler outer parts of 
the circumstellar disk (> 100 AU) and is spinning 
at a high enough rate. Thus, it is possible that 
some of our companions (and perhaps primaries) 
are unresolved binaries themselves. If that were 
the case, this could be reconciled with our common 
photometric distance estimates given the large as¬ 
sociated uncertainties and could be tested obser- 
vationally with either adaptive optics observations 
or high-precision radial-velocity monitoring of the 
primary. 


To differentiate the evolutionary states of our 
stars, we u se the four sp e ctral classes of metallicity 
defined bv iLmine et al.l ( 2007l l that are defined by 
the CTi 05 /CaH index, which is based on the CaH2, 
CaH3 and Ti05 molecular bands in the optical. 


This i ndex has been re-calibrated by IPhital et al 


(I2OI2II and is used to differentiate dwarfs, subd¬ 
warfs, extreme subdwarfs and ultrasubdwarfs. Us¬ 
ing this index for all of our M dwarfs, we found 
that they all have CriQb/CaH > 0.85, meaning that 
they all are in the dwarf metallicity class with 
near-solar metallicity. 


7. CONCLUSION 

We have discovered 14 new binary systems with 
companions of spectral types M6-L5 at separa¬ 
tions of 6"-170" from their primaries, correspond¬ 
ing to projected separations of 250-7500 AU at 
the distances of the systems. We also recovered 
nine already known binaries. Ten of our com¬ 
panions have a spectral type of LO or later, two 
of them being comfortably in the brown dwarf 
regime: 2M1115-bl607 (L5±l) and 2M1259-bl001 
(L4.5±0.5). The latter is a newly identified brown 
dwarf and orbits a mid-M dwarf. The most 
widely separated system is NLTT 687, consisting 
of an M3-I-L1 pair with a separation of 7400 AU. 
Other very wide systems are TYC 1725-344-1, a 
G5-I-M9 pair with a separation of 6700 AU and 


LSPMJ1441-I-1856, a M6-I-L1 pair with a separa¬ 
tion of 4110 AU. Pairs consisting of a G-type star 
and an ultracool dwarf (e.g., HD21746, G5-I-M6.5) 
provide an opportunity to cali brate the metallic¬ 
i ty sca le of M ultracool dwarfs Roias-Avala et all 
( 2 OIOII . This calibration is useful because deter¬ 
mining the metallicity of such M dwarfs would nor¬ 
mally require parallaxes and high-resolution spec¬ 
tra, which are expensive data to acquire and are 
limited to a few bright cool stars. While our sys¬ 
tems are not the most extremes, they can never¬ 
theless help better define the parameter space in 
which wide low-mass companions can form. 
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sity of Massachusetts and the Infrared Processing 
and Analysis Center/Galifornia Institute of Tech¬ 
nology, funded by the National Aeronautics and 
Space Administration and the National Science 
Foundation. This publication makes use of data 
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sion Laboratory/California Institute of Technol¬ 
ogy, funded by the National Aeronautics and 
Space Administration. This research has beneht- 
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Archive housed at http: //www. vlmbinaries. org 
and maintained by Nick Siegler, Chris Gelino, 
and Adam Burgasser. This research has bene- 
htted from the SpeX Prism Spectral Libraries, 
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ucsd. edu/~adain/browndwarf s/spexprism. 
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